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Abstract 

ON 

We calculated the polarization of ip' due to gg — ► QQ[ 3 fj 8 ^]<? — > color- 



octet quarkonia production. We find that at low transverse momenta the ij)' is 
unpolarized due to the contributions proportional to the L = S = and L = 



^ ■ S = 1 color-octet matrix elements. As p± increases, the ip' mesons become 

5h ; 

100% polarized, as predicted by fragmentation calculations. Polarization due 

to lowest order color-singlet production is also considered, which qualitatively 

has a similar shape to the color-octet production. 
14.40.Gx,13.87.Fh 
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The study of quarkonia production has recently received renewed attention due to new 
data and interesting theoretical developments. Traditionally, quarkonia production has been 
calculated in the color-singlet model (CSM), where the heavy quark- ant iquark pair is pro- 
duced in a color-singlet configuration at distance scales short compared to Aqcd- While 
the CSM is successful in describing many phenomenological aspects of quarkonia, it has 
become clear that it fails to provide a consistent picture of quarkonia production. Order of 
magnitude discrepancies have been found between CSM predictions and new measurements 
of ijj and T production at several colliders. These disagreements have called into question 
the validity of the CSM, and stimulated new ideas regarding quarkonia production. 

Quarkonia are inherently nonrelativistic due to the large mass Mq of the heavy quark and 
antiquark. Consequently, the physics of quarkonia involves a new small parameter, the ve- 
locity v of the heavy constituents inside the QQ bound state. An effective field theory called 
Nonrelativistic Quantum Chromodynamics (NRQCD) has been established based on a 
double power series expansion in the strong interaction fine structure constant a s = g 2 s /^ 
and the velocity v ~ 1/ logMg. NRQCD allows for the creation at short distances of a heavy 
quark- ant iquark pair in a color-octet configuration which later hadronized into a colorless 
final state quarkonia. Similar to Heavy Quark Effective Theory, NRQCD incorporates an 
approximate spin symmetry, which constrains various multiplet structures, transition rates, 
and polarizations. There have been many applications of NRQCD to quarkonia production 
in various high energy processes 0], but the validity of the picture still has to be verified. 

NRQCD makes definite predictions of the polarization of ^'s produced in a high energy 
collision. Thus, one test of the color-octet picture would be observing the polarization of ip 
mesons produced at the Tevatron consistent with NRQCD calculations. At large transverse 
momenta, quarkonia are primarily produced by gluon fragmentation The gluon is 

nearly real and transverse in the high p± limit, and the resulting QQf 3 ^ 8 ^] pair inherits 
this spin alignment. The long distance hadronization into a colorless ip preserves all angular 
momentum information, due to the NRQCD approximate spin symmetry. Thus, ip mesons 
produced at large p± are 100% transversely aligned [f|. Higher order a s correction to the 
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polarization of ip from gluon fragmentation have been calculated, and occur at the few 
percent level [[KJ. 

Gluon fragmentation is, however, only valid in the p± 3> Mq limit. At low transverse 
momentum, large numbers of ?/>'s are produced via color-octet states with L = S = and 
L = S = 1 [ [Tiyrm . Corrections to the fragmentation limit are not constrained to preserve 
the polarization of the if). Therefore, to use quarkonia polarization as a test of the color-octet 
mechanism, we need to investigate the spin alignment due to these states. In this paper, we 
will only consider the polarization of ip' mesons. A similar analysis for other charmonia and 
bottomonia states is possible, but these mesons are complicated by feeddown from higher 
level states. 

In the NRQCD formalism, the production cross section for a quarkonia state H in the 
reaction A + B — > H + X can be written as 

da(AB ^HX) = J2 C dx x dx 2 f a/A { Xl ) f b/B (x 2 ) £ da ab [n] (0%). (1) 

ab J ° n 

The first sum in Eq. ([]]) is over all partons in the colliding hadrons, and the parton distri- 
bution functions are denoted by f a /A and fb/B- The partonic cross section da^n] describe 
the production of a quark- ant iquark in a state n and can be calculated perturbatively in 
a s . The NRQCD matrix elements (O^) parameterizes the hadronization of state n into 
the quarkonium state H plus light hadrons [|l|. These NRQCD matrix elements contain 
all the nonperturbative information in the production process, and must be extracted ex- 
perimentally. The order in the velocity expansion at which each of these matrix elements 
participates in the ipQ creation processes is governed by NRQCD counting rules (T^j. 

The unpolarized cross sections for producing ip' mesons in quark- ant iquark, quark-gluon, 
and gluon-gluon scattering have been previously calculated up to 0(a^,v 7 ) [nj,F^]. At this 
order, the matrix elements that appear in Eq. (0) are (Of ( 3 *Si)), (0$ ( 3 Si)), (Of ( 1 5 , )), 
and (Of ( 3 Po))- The values of the color-octet matrix elements were extracted in JTjl^l by 
fitting the magnitudes of the calculated NRQCD cross section to Tevatron data, and are 

(Of CSi)) = (4.6 ± 1.0) x 10- 3 GeV 3 , (2a) 



{O *^ 0)) + { °*^ So)) = (5.9 ± 1.9) x KT 3 GeV 3 . (2b) 

Only the linear combination of (Of ( 1 >S'o)) and (Of ( 3 Po)) could be extracted. The error 
bars are statistical and do not reflect the large systematic uncertainties in heavy quark 
masses, color-singlet wavef unctions, parton distribution functions and next-to-leading order 



corrections ||14j| . These matrix elements should only be considered as an order of magnitude 
estimate, due to the large uncertainties. Color-singlet quarkonia production is suppressed 
relative to the color-octet production, as can be seen in Fig. 5 of Ref. [|rj . 



It is possible to use the results of Ref. [ 12 1 to obtain the polarization of the final state tp' 



when the intermediate QQ pair is in a 1 S , q 8 ' ) state. For the intermediate 3 S[ 8 ^ or 3 Pj 8 ' ) states, 
however, the analysis in Ref. 0] cannot be used to obtain a polarized cross section. As 



pointed out in Ref. [fTPf , since heavy quark spin symmetry is an approximate symmetry in 



the NRQCD Lagrangian, L z and S z are good quantum numbers. Therefore, to calculate the 
polarization of the final state ip' correctly, we must project the hard scattering amplitude 
onto states with definite L z and S z , square the amplitude, and then do the sum over L z . 
A formalism was developed later to correctly calculate the polarization in the NRQCD 
framework [15]. Also, since we are interested in the polarization of the quarkonium in the 
hadron frame, we cannot write the amplitudes completely in terms of partonic varibles.[] 
In Ref. the hard scattering amplitude was projected onto states with definite JJ Z and 



then squared, with the polarization vectors defined in the parton frame. This will give the 
correct unpolarized cross section |T6|JT5| , but will not be useful in calculating the polarized 



cross section. 



The methods used in calculating the amplitudes are similar to those described in Ref. [JT_2 
and the discussion will not be repeated here. The only difference being that we projected 
the amplitude onto states of definite L z and S z , squared, and then summed over L z . The 



1 We thank M. Beneke and M. Kramer for pointing out this error in the original version of this 
paper. 
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polarization vector e(A) of the quarkonium state then explicitly enters the differential cross 
section, Eq. (|I|), as fl4| 



-^M = A ab [n] + B ab [n](e(X) ■ h) 2 + C ab [n}(e(X) ■ k 2 f + D ab [n}(e(\) • fe 1 )(e(A) • h), (3) 

where k\ and k 2 are the momenta of the initial state partons a and b. The coefficients 
A, . . . ,D are shown in the Appendix. If we sum over the polarization vector, we recover the 



unpolarized cross sections from Ref. [12 



The ratio of longitudinal differential cross section to the unpolarized differential cross 
section, 

« = -rhr> (4) 

a T + a L 

can be measured in ip' — > £ + £~ decay. The leptons are distributed in angle according to 



dcos9 



oc 1 + a cos 6, (5) 



where 



l-3£ 

a = Y+^' (6) 

and ^ denotes the angle between the lepton momentum in the ip' rest frame and the ip' 
momentum in the lab frame. In Fig. 1, we plot a for prompt ip' production at the Tevatron.0 
Since there is only a value for the linear combination in Eq. (f2~b|) , we cannot give a definite 
prediction for a. Instead, the solid curve represents a when (Of ( 3 -Po)) = 0. The dashed 
curve illustrates a when the contribution from (Of ( 1 5'o)) is set to zero. The shaded region 
illustrates the effect of the uncertainties in the matrix elements in Eq. (0). 

The angular distribution approaches the transverse form 1+cos 2 9 at highp^ as predicted 
by gluon fragmentation computations II . At low transverse momentum, the ip' is essentially 



2 In this paper, MRSDO parton distribution functions evaluated at the renormalization scale \i = 
/p j_ + 4M| were used, with M c = 1.48 GeV. A pseudorapidity cut of \rj\ < 0.6 was imposed. 
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FIG. 1. Coefficient a which governs the lepton angular distribution in ip' — » £ + £~ decay plotted 
as a function of p± . The solid and dashed curves illustrate a for ip' production at the Tevatron 
when (Og ( 3 Po)) and (0§ ( 1 5'o)) respectively vanish. The shaded region shows the effect of the 
uncertainty in the extraction of the matrix elements. The dotted line corresponds to the lowest 
order, color-singlet production. 

unpolarized due to L = S = and L — S — 1 color-octet states. Since the two curves are 
similar in shape, the true value for the angular coefficient should be close to the curves 
shown. The effect of uncertainties in the matrix elements and higher order corrections can 
be qualitatively described by slight displacements of the curves in Fig. 1, without changing 
the asymptotic behaviors. Higher order a s corrections, however, can change the asymptotic 
high p± behavior, but it should be a small effect [[U] . 



Also plotted in Fig. 1 is the polarization due to the lowest order, color-singlet production, 
gg — > ccpsj 1 ^. As can be seen, the qualitative shape of the color-singlet curve is similar to 
the color-octet curves. While higher order color-singlet corrections will modify this shape, 
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it is clear that by just observing the qualitative shape presented in Fig. 1 of the prompt ip' 
polarization at the Tevatron will not be a clear signature for color-octet production. 
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APPENDIX: 



Functions entering into polarized cross sections: 3 
qq — > ip^g: 

4a 3 vr 2 (4(t 2 + m 2 ) - iuj (s 2 - 2tu + M A ) 



3 e( 8 )i 



Aqq [ S 1 



81M 3 s 2 tu(s - M 2 ) 2 



R psf) = _i!^ZL I! ^ / (Ala) 

" L 1 J 81Ms 2 tu(s-M 2 ) 2 ' 1 ; 



3 0(8)1 



16a 3 vr 2 (4(t 2 + u 2 ) - iu 



" [ 1 J 81Ms 2 tu(s-M 2 ) 2 
D qq [ s S{ 8) }=0, 



gq — > ip^q: 



3 (8) _ 80a 3 7r 2 g 2 -2t£ + 3M 



27M 3 s 2 


s(s-M 2 ) 2 ' 


640a 3 7r 2 


£« + uM 2 - M 4 


27Ms 2 


s 2 (s - M 2 ) 3 


640a 3 7r 2 


£u + tM 2 - M 4 


27Ms 2 


s 2 (s-M 2 ) 3 ' 


640a 3 7r 2 s 


2 + sM 2 - 2tu 



R r3 P (8)i_ Q4ua^ m + w-M- 



qq[ J J 27MS 2 s 2 {s-M 2 f ■ 



4 r3 c(8) 1 _ <*W (4(g 2 + U 2 ) ~ SU) (g - 2g£ + M 4 ) 

99 [ 1 J 54M^ 2 sw(t-M 2 ) 2 

3 (8) 2a 3 7r 2 (4(g 2 + u 2 ) - su) 

^ J " 27MI 2 " su(t-M 2 ) 2 ' (A2a) 

r r3 0(8)1 4ag7r 2 (4(g 2 + £ 2 ) - su) 

Dl J 27Ms 2 <m(t-M 2 ) 2 ' 

n [80(8)1 _ 4 ^ 3?r2 (4(g 2 + m 2 ) - g^) 

ot[ 1 J 27M5 2 sn(t-M 2 ) 2 ' 



3 The coefficients corresponding to (0g ( 1 5'o)) production can be obtained from Ref. \\2\ and will 
not be repeated here. 
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99 



A 



ii'i 



B 



>i<i 



a 



91 



D 



.1 



M 8, l 



3p(8)! 



</'/ 



3p(8)i 



10aV i 2 - 2su + 3M 4 



3p(8)i 



9M 3 s 2 t(t - M 2 ) 2 


80a 3 vr 2 


sit + itM 2 - M 4 


9Ms 2 


t 2 (t-M 2 ) 3 ' 


80a 3 vr 2 


t + M 2 


9Ms 2 


t 2 (t-M 2 ) 2 ' 


80a 3 vr 2 


£ 2 -M 2 (2s + £) 



9Ms 2 t 2 {t-M 2 f 



Agg 


[ 3 sf] 


10a 3 s 7i 2 M 


s 2 {s - 


M 2 ) 2 + (si + iu^ 


- MS 


-s 2 ) 


8 


Is 2 


(s 


-M 2 ) 2 (t-M 2 ) 2 (u- 


•M 2 ) 2 


B 99 


[ 3 M 1} ] 


20a 3 7r 2 M 3 


(s 2 + t 2 ) 








81s 2 


(s- 


- M 2 ) 2 (i - M 2 ) 2 (u - 


M 2 ) 2 ' 


Cgg 




20a 3 7r 2 M 3 


{s 2 + u 2 ) 








81s 2 


(s- 


- M 2 ) 2 (t - M 2 ) 2 {u - 


M 2 ) 2 ' 


Dgg 


[ 3 M 1} ] 


40a 3 vr 2 M 3 


s 2 








81s 2 


(§- 


- M 2 ) 2 (t - M 2 ) 2 (ii - 


M 2 ) 2 ' 


A 99 


[ 3 ^ 8) ] 


n ? 


7T 2 


27(s 2 


- iu - M 2 s) + 19M 4 






36M 3 s 2 ( 


s - M 2 ) 2 (t - M 2 ) 2 (u - M 2 ) 2 








X 


s 2 (s- 


-M 2 ) 2 


+ tii (si + iii + us — 


s 2 )_ 




B 99 


[ 3 sj 8) ] 


a 3 7r 2 


27(s 2 


- iu - M 2 s) + 19M 4 


(s 2 + t 2 ) 


U 




is- 


- M 2 ) 2 (i - M 2 ) 2 (ii - 


M 2 ) 2 


Cgg 


3 5 (8)j 


a 3 7r 2 


27(s 2 


- iu - M 2 s) + 19M 4 


(s 2 + it 2 ) 




,Ms 2 




- M 2 ) 2 (i- M 2 ) 2 (ii- 


M 2 ) 2 


Dgg 


[^i 8) ] 


a 


> 2 " 


27(s 2 - 


- iu - M 2 s) + 19M 4 


s 2 





(A2b) 



(A3a) 



(A3b) 



We have introduced the variable z = V tit to simplify some of the coefficients. 
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A gg [ 3 Pf ] ] = ^|^|m 2 s 3 (s - M 2 ) 3 (s 4 - 2M 2 s 3 + 7M 4 s 2 - 6M 6 s + 3M 8 ) 

+s 2 z 2 (s - M 2 )(s 6 - 8M 2 s 5 + 23M 4 s 4 - 50M 6 s 3 + 56M 8 s 2 

-31M 10 s + 6M 12 ) 
-H 4 (4s 6 - 9M 2 s 5 + 31M 4 s 4 - 71M 6 s 3 + 77M 8 s 2 - 34M 10 s + 6M 12 ) 
+5 6 (6s 5 + 4M 2 s 4 + 20M 4 s 3 - 33M 6 s 2 + 22M 8 s - 3M 10 ) 
-2i 8 (2s 3 + 2M 2 s 2 + 5M 4 s - 3M 6 ) 
+5 10 (s - M 2 )}/ (sz 2 (s - M 2 f{M 2 s + 5 2 ) 3 ) , 

^[ 3 Pf] = -^{^ 5 (M 2 -^ 

-iu 3 (M 2 - m) 4 (M 8 - 7M 6 « + 42M 4 m 2 - 52M 2 m 3 + 24m 4 ) 

+i 2 u 2 (M 2 - m) 3 (2M 10 - M 8 u - 39M 6 m 2 + 152M 4 m 3 - 166M 2 m 4 

+68m 5 ) 

-t 3 u(M 2 - m) 2 (M 12 + 9M 10 « + 2M 8 m 2 - 134M 6 m 3 + 361M 4 m 4 

-339M 2 m 5 + 116m 6 ) 
+i 4 u(M 2 - m)(11M 12 + 9M w u + 16M 8 m 2 - 274M 6 m 3 + 589M 4 m 4 

-471M 2 w 5 + 128m 6 ) 
+t 5 (M 2 - w)(4M 12 - 51M 10 w + 2M 8 w 2 - 36M 6 w 3 + 282M 4 m 4 

-329M 2 m 5 + 80m 6 ) 
-t 6 (20M 12 - 129M 10 m + 94M 8 m 2 - 12M 6 w 3 + 150M 4 m 4 - 147M 2 m 5 
+8m 6 ) 

+8f (5M 10 - 19M 8 m + 6M 6 m 2 + 6M 4 m 3 - 3M 2 m 4 + 5m 5 ) 
-8t 8 (5M 8 - 11M 6 m - 2M 4 m 2 + 7M 2 m 3 - 5m 4 ) 
+20t 9 (M 2 -m) 2 (M 2 + m) 

-4t 10 (M 4 - m 2 )}/ (sH 2 u 2 (s - M 2 f(t - M 2 f(u - M 2 f) , 
C gg ?P ( P]=B 9g fpf\^ (A3c) 
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-M 2 s A z 2 (s - M 2 ) 2 (22s 3 - 38M 2 s 2 + 19M 4 s - 4M 6 ) 
-2s 3 z 4 (s 5 - 22M 2 s 4 + 62M 4 s 3 - 62M 6 s 2 + 27M 8 s - AM 10 ) 
+s 2 z 6 (2s 4 - 17M 2 s 3 + 66M 4 s 2 - 31M 6 s + 8M 8 ) 
+2H 8 (3s 3 - 6M 2 s 2 - 3M 4 s + 2M 6 ) 

-2^ 10 (5s - 3M 2 ) + 4S 12 } / (s 2 ^ 4 (s - M 2 f(M 2 s + 5 2 ) 3 ) . 
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